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Abstract cavity wall, around the i\iel and is then exhausted 

through the nozzle. The objective of this study 


The uranium plasma nuclear rocket is a pro- 
posed device which features a high specific im- 
pulse. The objective of this study was to find the 
upper limit of the ability of the seeded hydrogen 
propellant to intercept the thermal radiation from 
the high temperature fuel and thus protect the 
reactor cavity wall from excessive heat flux. This 
limit fixes the maximum reactor power that can be 
used for a fixed propellant flow rate; and deter- 
mines the maximum achievable cavity specific im- 
pulse. This maximum cavity propellant temperature 
can then. be used to find the amount of nozzle wall 
seeded transpiration flow that is necessaiy. This 
required nozzle coolant flow reduces the maximum 
cavity specific impulse to some lower value that is 
the maximum specific impulse such as an engine sys- 
tem can provide. The model used for the heat- 
transfer calculation is as follows: flow enters 

the cavity wall at a given temperature, and picks 
up the amount of radiant heat deposited on the 
inner surface of the wall as it flows through it. 
The propellant then leaves the wall at the solid- 
wall -limited temperatxire and begins to flow radi- 
ally inward inside the cavity. Between the wall 
and the central uranium plasma there is an energy 
balance between the absorbing hydrogen flow toward 
the plasma and the radiant heat flux toward the 
wall. The same basic heat-transfer analysis is. 
used in the cavity and in the nozzle. However, in 
the nozzle, one must additionally specify the 
boundary layer thickness caused by the injection of 
the absorbing, transpiration coolant flow. The 
thickness of this coolant layer in the nozzle is 
calcilated from conservation of mass and momentum. 

This study showed that, for an eight foot 
cavity diameter uranium plasma nuclear rocket op- 
erating at a pressure of 1000 atm and a propellant 
mass flow rate of 10 Ibm/s, the cavity wall could 
be cooled up to a power level of 7400 MW. This 
corresponded to a maximum cavity specific impulse 
of 5800 s. The wall heat flux was very low for 
any reactor power below this limit of 7400 MW. The 
reason for this was that there was a relatively 
cool, opaque insulating layer of seeded propellant 
between the hot plasma and the solid wall. The 
degradation effect of the transpirationally cooled 
nozzle on the cavity specific impulse resulted in 
a maximum specific impulse of 5200 seconds at a 
reactor power of 7400 MW. Approximately 12 percent 
of the total propellant was used to cool the 
nozzle wall. 

Introduction 

The uranium plasma nuclear rocket is a pro- 
posed device which features a high specific im- 
pulse (2000-7000 sec) and a moderately high thrust 
(10,000-100,000 Ib).(l) A conceptual design of 
current interest is shown in Fig, 1 and is called 
a "spherical" or "porous-wall" reactor. Energy is 
generated by nuclear fission in a uranium plasma 
fuel. This energy is radiated to a seeded hydro- 
gen propellant. The propellant flows through the 


was to determine the upper limits on specific im- 
pulse placed by the wall cooling requirements of 
both the reactor cavity and the exhaust no 2 ^ 1 e. 

As the propellant flows through the cavity 
wall, it is heated from cryogenic temperatures to 
solid-wall -limited temperatures. Thus, the pro- 
pellant transpirationally cools the cavity wall. 
The amo\mt of energy deposited in the propellant 
as it flows through the wall equals the amount of 
energy deposited on the cavity side of the wall by 
the impinging radiant eind conductive heat flux. 
(This assumes the gamma heating of the wall is 
small. ) Thus, specifying inside and outside wall 
temperatures (and the mass flow through the cavity 
wall) fixes the allowable cavity wall heat flux. 
Using the radial component of the energy equation, 
a modified diffusion approximation for the radiant 
heat flux and an assumed radial mass flow profile 
inside the cavity, the radial temperature profile 
and heat flux profile from the edge of the fuel to 
the cavity wall can be calculated. From these 
profiles the reactor power, the cavity Isp, and 
the thrust can be calculated. 

As the propellant leaves the reactor cavity 
and flows throu^ the nozzle, it begins to heat 
the nozzle wall by forced convection and by therm- 
al radiation. A layer of seeded coolant gas must 
be injected through the nozzle wall. This coolant 
transpirationally cools the nozzle wall in much 
the same way as the propellant cools the cavity 
wall. In the nozzle the he at -transfer dimension 
is much smaller than in the cavity; therefore, 
the coolant flow rate per unit area must be much 
larger in the nozzle than in the cavity. 

The approach used in this stuc3y was first to 
find the cavity wall cooling limitation; this will 
determine the maximum reactor power and the maxi- 
mum cavity specific impulse. Then, using the cor- 
responding propellant outlet temperature, the 
amount of required nozzle wall coolant is deter- 
mined, and the cavity specific impulse is reduced 
due to this additional flow of cooler gas. The 
resulting specific impxilse is the maximum that can 
be produced by the engine without burning out 
either the reactor cavity wall or the nozzle wall. 

Symbols 

A area, m^ 

a absorption coefficient, m“^ 


En expodential integral of order n, Ej^(z) = 



h enthalpy, J/kg 
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I ■; intensity^ w/Cm^-sr-Hz) 

t-'hermal.. conduct ivity, w/m 
,L ■ cooling surface length, m 
M ■■■ ■ mass flow rate, kg/s 
.n '■ normal ;unit vector 
. n ■■■ ■: index \ised: in Taylor Series 
P' i ■ 

q[- heat' flux,- w/m^ 

: r VI i radius,; m 

Vt V ;-. ^&r^eratuj’e,- K' [ ' ' ' 

':v ';l;r velocity, m/s ■ • 

Vz. i'^'-'V ^^iaJ“/.length, mV 

5 . distance from solid wall towards plasma, m 
\0 ' nozzle,; rad 

:v ■' VvV?^^l^®^^^> . Hz ' . '* ' ■ 

; p .kg/m^ 

Q : .Stefen-Boltzmann constant, 5.6697x10”® w/ 

T ..optical, depth from solid wall towards plasma 

0 ^ . ■ r^olid au.gle , .sr 

VS.x]hs,efiptsr 

:,c Vv V cpnduet.ipnV ; 

■ ■•V'Vhp^zle" - 

o . cool side Qf-.solidwall - 

;'p-VV';-'pkasma--; 

;R-;il V.fadiatxon ;■ 

Vt; ■:;:V;hoz:zle''thrpat ■ 

plasma- side;; of' wall- ; 

;0--' :;i;Vup;Vstream;- of -nozzle throat- ^ ^ 

Analysis 

. . The model, used for the heat-transfer calcu- 
lation is shorn in Pig. 2. Plow enters a wall at 
a temperat'ui’e, T^, then picks up the amount of 
heat deposited on the wall, and then leaves the 
waa,l at a temperature, T^. In terms of enthalpy, 
:.h,i^this:-; heafihaian 

= -P-^(1V “^ ^o) (1) 

; VInside th^^ there is an energy /balance 


between the radially inward convection toward the 
plasma and the radially outward heat flux toward the 
wall. In vector notation this can be written as; 

0 = (^) • Vh + V . ^ (using only the radial terms) 

( 2 ) 

The heat flux can be written as the sum of two 
terms; conduction and radiation, 

<1 = <lc + (3) 

The conductive heat flux is proportional to the 
local temperature gradient, 1 

/ ^ (4) 

The radiative heat flux is given in Ref, 2 as 
"^R V= ^ n 

For a grey gas in local thermodynamic eg.uilibr iron, 
in cartesian coordinates, the radiative heat flux 

is .". 

% - - 2aT^3(Xp- x) ; 

+ r"^ 2aTtE3(x - x')dx' - /" V' 2aT%2('T' - x)dT' 

/: *0.. 

( 5 ) 

An approximation to this heat flux^ can be made if 
the radiation from the propellant can be neglected. 
Thus'--'if 

2aT|E3('^) - al^» ctT^( 1 - 2E3(xp)) 

then ■■ ■ 

U = 20T^3(t) - 2oT^3(Xp - x) (5a) 

(This approximation is used when is small.) 

Another; approximation to this heat flux relation 
can be made if the source function; for ^ the gas can 
be represented by a Taylor Series , Thus 

T^= T-^x) -t- — (x' - x) + . . . : 

dr 

(x*; - x)” i/t/t) 

■■ : n! • ■ : an:’^ 

n=o : ^ 

Substituting this into Eq. (5); and; integrating 
gives; 
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- 2orjE3(Tp - t) 

n=0 U=0 

1 . (1 - (-1)^ 
J (n+2) 


Tp-r) 




- E (t)) + 2 c ^ 

d.T 


[' 


(Tp - T)E3(Tp - t) 


Combining Eqs. (7) and (8) yields 

6 2 (pv)^ L/yp AP* 

An analogous equation for a convergent -divergent 
nozzle, Fig. 4, is 


Neglecting the terms of the series for which n is 
greater than one yields 

qp = 2aT^3(T) - 2aT^3(Tp - t) + Ecrl'^CEjCTp - t) , 


( 9 ) 


5 ^ 




where Eqs. (8) and (9) have been rewritten in cyl- 
indrical geometry, and the control volume is a 
convergent -divergent shell next to the nozzle wall. 
The total mass flow through the nozzle wall is 


+ tE 3 (t) + E 4 (Tp " t) + E 4 (t) - 2/3J 

If the optical thickness (optical depth from the 
wall to the radiating plasma) is large. 


Tp »1 

then 


2a FtEsCt) 
dT L 



If the local propellant optical depth is large. 


then 


T » 1 


qR = - 


3 dT 


(5c) 


o 

% = 


^■\)/\ ~ ^ (AyA^ - 1) 

sin 6q sin Sq 


(P'^)w H (^) 


Properties for uranium vapor were taken from 
Refs. 3 and 4. Properties for hydrogen were taken 
from Refs. 5, 6, 7, and 8. The absorption co- 
efficient for the solid, seed (uranim 238) was 
assumed to be 50,000 cnr/g. 

Discussion 


A cross section of the uranium plasma nuclear 
rocket is shown on Fig. 1. The propellant flows 
through the cavity wall, transpirationally cooling 
the wall. Then the propellant flows from the wall 
toward the fissioning uranium plasma convecting the 
radiant heat that it absorbs back toward the fuel. 
Near the fuel the propellant turns and flows 
tangentially out of the cavity by way of the nozzle. 
In the nozzle the propellant acts as a source of 
radiant heat. Additional coolant must be used to 
transpirationally cool the nozzle wall and absorb 
the radiant heat. 


which is just the standard diffusion approximation. 

The optical depth is related to the distance 
from the cavity wall by 

dT = a d5 (6) 

Using the given wall temperature, the calculated 
wall heat flux from Eq. (l), and the given wall 
distance as boundary conditions, the differential 
Eqs. (2), (3), and (6) coupled with the auxiliary 
Eqs. (4), (5) and an assumed variation of mass 
flow rate, (pv(b) equals a constant), can be inte- 
grated numerically from the wall into the cavity. 

The heat-transfer analysis used in the nozzle 
is the same as that used in the cavity. A flow 
model used in the calculation of the coolant 
boundary layer thickness is shown in Fig. 3. Con- 
servation of mass requires 

(pv)^ L = (pv)3_ 5 (7) 

Conservation of momentum, (neglecting compressible 
and viscous effects) requires 

AP S ^ (pv^)l 6 (8) 


Figure 5 illustrates the complexity of the 
problem. It shows the "Rosseland Mean Absorption 
Coefficient" versus temperature. At low tempera- 
tures the solid seed is the dominant absorber. As 
the temperature increases, the solid- seed absorp- 
tion coefficient decreases because the propellant 
density decreases and, therefore, the number of 
solid particles per unit volume also decreases. 

At some temperature (~5000° K) the solid seed be- 
gins to vaporize and the absorption coefficient de- 
creases very sharply. In the temperature region 
between 7000-10,000° K, the absorption is due main- 
ly to seed vapor which in this case is uranium 
vapor, since in these calculations we have used 
depleted uranium as the seed material. At temper- 
at-ures above 10,000° K, the absorption is due main- 
ly to the hydrogen becoming absorptive. 

Figure 6 shows results of some typical reac- 
tor cavity calculations. Figure 6(a) shows the 
Isp and inside cavity wall temper atiure versus 
reactor power for a fixed propellant flow rate of 
10 Ibm/s, cavity diameter of 8 ft, and seed concen- 
tration of 10 percent by weight. The Isp is 
proportional to the square root of the reactor 
power. For a reactor power below about 7400 MW, 
the cavity inside wall temperature remains nearly 
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constant at 944^ K. At a reactor power equal to 
about 7400 the wall temperature has increased 
rapidly to a value that approaches and would soon 
exceed solid wall temperature limits. This value 
of reactor power would be the maximum allowable 
reactor power, and for this particular set of input 
parameter values corresponds to a maximum attain- 
able cavity specific impulse of 5800 seconds. Be- 
cause reactor power could fluctuate, a "nominal 
operating power" level would probably be fixed at 
90 percent or so of the maximum operating power. 

Figure 6(b) shows a radial temperature profile 
for both the maximum and the nominal: reactor power 
situations. Both temperature profiles are similar 
except near the cavity wall. For the nominal case 
there is a cold opaque layer of propellant between 
the . hot propellant and the solid wall. For the 
maximum, ease: this layer has been; reduced to almost 
zero thickness. ' 

The model: used to estimate the wall cooling 
iz^quirementsiof the nozzle is shown in Fig. 4, 
Bi>opellant flows from the cavity (which is at a 
pre.ssure of 1000 atm) through the nozzle, and is 
discharged into space. Transpiration cooling with 
an absorbing gas is necessary in the high acceler- 
ation (throat) region of the nozzle, a region 
taken to correspond to local static pressure vary- 
ing from 900 atm to 100 atm. The thickness of the 
injected boundary layer gas was calculated from a 
momentum balance flow model assuming negligible 
viscous shear. 

Figui'e 7 shows the result of nozzle cooling 
limitations on Igp. It shows .the nozzle Igp 
(calculated using the propellant flow from cavity 
plus the nozzle coolant flow) versus cavity Igp 
for a seed concentration of 40 percent and a wall 
temperature of 1100^^ K.; The spread in the results 
is due to uncertainties related to what extent the 
nozzle coolant would mix down stream of the nozzle 
throat with the main propellant.: The cavity-wall- 

limited Igp for this case is about 5800 seconds, 
while the actual Igp including nozzle cooling 
degra-dation is between 5100 and 5500 seconds. 

These results are not intended to be the final 
word on the complicated subject of Igp limita- 
tions of a gas -core rocket engine due to thermal 
protection requirements. These results were ob- 
tained by a "first-cut" look at the problem. The 
answers appear encouraging, but more information 
on seeded hydrogen opacity and on transpiration 
flow- in the presence of themal radiation absorp- 
tion is required. 

Conclusions 

shown that for an eight foot 
cavity diameter uranium plasma nuclear rocket op- 
erating at a pre ssure of 1000 atm, and a propel- ; 
lant mass flow rate of 10 Ibm/sec, the cavity wall 
could be cooled up to a reactor power level of 
7400 M. This corresponded to a maximum cavity 
specific impulse of 5 s. The wall heat flux 
was very low: for any reactor power below this limit 
of 7400 MW. The reason for this was that there 
was a relatively cool, opaque insulating layer of 
: seeded propellant between the hot plasma and the 
solid wall. 

The additional coolant flow required to pro- 


tect the nozzle reduced the maximum cavity specific 
impulse value of 5800 seconds to a value of 5200 
seconds. Approximately 12 percent of the total 
propelXant was used to cool the nozzle wall. 
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Heat transfer model of spherical reactor. 
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Figure 3 





Figure 5. - Rosseland absorption coefficient 
. versus temperature at 1000 atmosphere, 

1. 01x10'^ N/m^. 


TEMPERATURE, K SPECIFIC IMPULSE, 



Reactor power, meg-W 

(A) SPECIFIC IMPULSE AND WALL TEMPERATURE 
VERSUS REACTOR POWER FOR A PRESSURE OF 
1000 ATM, LOlxlO+S N/m^; A PROPELLANT 
FLOW RATE OF 10 LBM/S, 4.54 KG/S; A CAVITY 
RADIUS OF 4 FOOT. 1.22 m; A FUEL/ CAVITY 
VOLUME RATIO 0.3; AND A MASS FRACTION 
OF SEED 0. 1. PROPELLANT ENTERS THE WALL 
AT 944 K. 



(B) CAVITY PLASMA TEMPERATURE VERSUS RADIUS, 
AT 1000 ATMOSPHERES, 1. 01x10+8 N/m2, FOR 
MAXIMUM AND NOMINAL REACTOR POWER; PRO- 
PELLANT-FLOW RATE =10 LBM/S, 4. 54 kg/S. 
PROPELLANT ENTERS THE WALL AT 944 K. 


Figure 6 


WALL TEMPERATURE, 
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SPECIFIC IMPULSE OF CAVITY, S 

Figure 7. - Specific impulse degradation due to 

nozzle coolant flow, (40% + 60% H 2 ). 

Coolant enters the wall at 300 K and leaves the 
wall at 1100°K. 
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